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Abstract 

Disulfide-rich peptide toxins found in the secretions of venomous organisms such as snakes, spiders, scorpions, leeches, and 
marine snails are highly efficient and effective tools for novel therapeutic drug development. Venom peptide toxins have 
been used extensively to characterize ion channels in the nervous system and platelet aggregation in haemostatic systems. 
A significant hurdle in characterizing disulfide-rich peptide toxins from venomous animals is obtaining significant quantities 
needed for sequence and structural analyses. Presented here is a strategy for the structural characterization of venom 
peptide toxins from sample limited (4 ng) specimens via direct mass spectrometry sequencing, chemical synthesis and NIVIR 
structure elucidation. Using this integrated approach, venom peptide Tvl from Terebra variegata was discovered. Tvl 
displays a unique fold not witnessed in prior snail neuropeptides. The novel structural features found for Tvl suggest that 
the terebrid pool of peptide toxins may target different neuronal agents with varying specificities compared to previously 
characterized snail neuropeptides. 
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introduction 

There is a need to discover new compounds that enhance the 
drug pipeline for treating ailments. The disulfide-rich peptide 
toxins found in venomous organisms such as snakes, spiders, 
scorpions, leeches, and marine snails are highly efficient and 
efTective for manipulating physiological pathways and for novel 
therapeutic and insecticide drug development [1-4] . In an effort to 
identify novel disulfide compounds from venomous snails, we have 
investigated the Terebridae, a sister group to the more familiar 
cone snails. Novel peptide toxins from the Terebridae would 
enhance the pool of therapeutic compounds available for 
biomedical research. Venomous disulfide-rich compounds include 
ion channel and enzyme inhibitors, growth factors, and structural 
or ligand-binding domains [5]. The disulfide-rich peptide library 
produced by venomous organisms is estimated to contain well over 
2 million compounds. The potential and promise of these 
compounds are evidenced by breakthroughs such as, angioten- 
sin-converting enzyme (ACE)-inhibitor, CaptoprU [6] and more 
recently, ziconotide (Prialt), the first drug from a venomous marine 
snaU, Conus magus, which is used to treat chronic pain in HIV and 
cancer patients [7]. Peptides, while currently not as bio-available 



as small molecule therapeutics, have advantages of higher target 
specificity and selectivity, and decreased toxicity. Despite then- 
advantages, the use of peptides in drug development lags in the 
discovery and characterization stages. Obtaining workable quan- 
tities of venomous natural products has remained a significant 
challenge due to the large amount of material that has traditionally 
been needed to identify and structurally characterize peptidic 
compounds [8-12]. A sensitive and robust method for discovery 
and characterization of disulfide-rich peptides from sample limited 
venom sources would facilitate their application as biochemical 
tools and potential drug development targets [13-15]. Presented 
here is an integrated approach of liquid chromatography mass 
spectrometric sequencing, chemical peptide synthesis, and NMR 
for the rapid identification and characterization of sample limited 
bioactive disulfide-rich venomous snail peptide toxins (Fig. 1). 
Structural characterization of Terebra variegata peptide toxin Tvl 
was accomplished using a total venom sample of 4 ng, an amovmt 
that would previously make characterization of venom compo- 
nents almost impossible [16,17]. 

Tv 1 is the first peptide structurally characterized from a terebrid 
snail. Terebrids are part of the Conoidean superfaniUy of 
predatory moUusks, which includes cone snails and turrids 
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Figure 1 . Discovery and cKiaracterization of disulfide-rich Tvl teretoxin. An integrated approach for characterizing sample limited disulfide 
peptidic natural products was applied: a) RP-HPLC on-line separation of Terebra variegata venom (highlighted is the region where the sequenced 
peptide eluted). b) ETD and CAD MS/MS analysis recorded on the native peptide to determine the sequence of Tvl peptide, c) RP-HPLC of linear and 
folded versions of chemically synthesized Tvl peptide, d) MS/MS spectrum of partially reduced and alkylated Tvl used to determine the disulfide 
connectivity, e) NMR solution structure of folded Tvl peptide, f) Bioactivity of 20 \iM of Tvl (-) versus normal saline (NS) (— ) solution injected into 
polychaete worms. Y-axis indicates level of activity, where 3 is the normal activity of the NS injected worm, 1 and 2 indicate decreased activity due to 
partial paralysis. 

doi:10.1371/journal.pone.0094122.g001 



[18,19]. Tvl is a novel twenty-one amino acid teretoxin peptide 
with a cysteine scaffold similar to the M-superfamHy of cone 
snaU neurotoxins, CC-C-C-CC. Apart from the cysteine 
pattern, there is no sequence homology between Tvl, TR(I/ 
L)CCGCYWNGSKDVCSQSCC, and known M superfamily 
conotoxins from Corns snails (Table 1). Conoidean venoms are a 
rich source of neuroactive disulfide-rich peptides, used to 
investigate cellular communication in the nervous system [17- 
22] and have diverse biomedical applications, including pain 
therapy [23,24]. In contrast to conotoxins, terebrid toxins, 
teretoxins, which were recently identified, [25,26] are an untapped 
resource for discovering novel neuropeptides. The sequence and 
structural characterization approach presented here significantly 
reduces the number of snails needed, while allowing for rapid 
characterization of venom components, in an effort to conserve 
natural resources (Fig. 1). 

Results and Discussion 

Venom Extraction and De-Novo Sequencing of Tvl 

The primary sequence characterization of Tvl was accom- 
plished using charge-enhanced Electron Transfer Dissociation 
(ETD) in combination with traditional CoUisional Activated 
Dissociation (CAD) on reduced and alkylated Terebra variegata 
venom [27]. The peptide toxins in the crude venom were reduced 
and subsequently alkylated using either iodoacetamide or N, N- 
dimethyl-2-chloro-ethylamine [27]. The latter converts every 
cysteine residue into a dimethyl lysine analogue, which increases 
the charge state of each cysteine-rich peptide toxin, resulting in 
excellent sequence coverage with ETD. 

Combining complementary dissociation strategies CAD and 
ETD facilitates de novo sequencing, as each dissociation strategy has 
specific fragmentation characteristics, which can be used to 



confirm the corresponding ETD and CAD spectra (Fig. 1). This 
is advantageous as peptide toxin sequences can be obtained 
directly from crude venom without prior fractionation and using 
only mass spectrometry [2 7] . This feature is particularly important 
for sequencing peptide toxins from animals producing minimal 
amounts of venom, such as terebrids, where only nanogram 
quantities are extracted from each snail specimen. Four nano- 
grams of crude Terebra variegata salivary gland contents were 
extracted with 30% acetonitrile/water in 0. 1 % trifluoroacetic acid. 
An aliquot of the crude venom was analyzed using LC-MS on an 
Orbitrap XL mass spectrometer to determine the number of 
teretoxins present. LC-MS analysis was repeated after reduction of 
the crude venom to determine the number of cysteine residues 
present in each component. At this point, Tvl was selected for 
further analysis due to its similarity in molecular weight and 
number of cysteines to M superfamily conotoxins. A reduced 
version of Tvl was dissociated using CAD (Fig. lb, and Figure SI 
in File SI) and the dimethyl lysine analogue version of the 
teretoxin was dissociated using ETD (Fig. lb and Figure S2 in File 
SI). The ETD MS/MS spectrum resulted in almost 100% 
sequence coverage. The sequence assignment by ETD could be 
supported with the CAD spectrum on the reduced toxin (Fig. lb 
and Figxire SI in File SI). There were two places of ambiguity in 
the Tvl peptide sequence: (1) the inherent uncertainty of 
isoleucine (He) and leucine (Leu) at the third position, which 
cannot be differentiated with the mass spectrometer used for this 
analysis, and (2) the order of the first two amino acids, Thr, Arg or 
Arg, Thr. The second ambiguity was resolved based on the 
absence of the cl ion, and the resulting Tvl sequence was 
determined as TR(I/L)CCGCYWNGSKDVCSQSCC. However, 
the He /Leu ambiguity at position 3 could not be resolved as there 
was not sufficient crude extract remaining after MS characteriza- 
tion. 
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Table 1. Comparison of Teretoxin Tvl and known M superfamily Conotoxins. 





Teretoxin 


Sequence 


Target(s) 


Reference 


Tvl 


TRICCGCYWNGSKDVC5QSCC 


unknown 


This work 


Conotoxin 


Mrllle 


VCCPFGGCHELCYCCD* 


unknown 


[23,33,49] 


Txllla 


CCSWDVCDHPSCTCCG* 


unknown 


[50,51] 


Mrllla 


GCCGSFACRFGCVOCCV 


unknown 


[50] 


Txlllb 


CCPPVACNMGCKPCC * 


unknown 


[31,50] 


Reg 12a 


GCCOOQWCGODCTSOCC 


unknown 


[52] 


Tx3f 


RCCKFPCPDSCRYLCC* 


unknown 


[31] 


Qc3.1 


ACCDPDWCDAGCYDGCC 


unknown 


[23] 


BulllA 


VTDRCCKGKRECGRWCRDHSRCC* 


Na^l .2,1.3, 1.4 


[53] 


BulllB 


VGERCCKNGKRGCGRWCRDHSRCC* 


Na^l. 2,1. 3,1 .4,1. 5,1. 7 


[53] 


BulllC 


IVDRCCNKGNGKRGCSRWCRDHSRCC** 


Na„1.2,1.3,1.4,1.7 


[53] 


MIMA 


ZuCCNVrNCjCbtjKWCKUHAUCC* 


TTX-r 


[54] 


KIIIA 


CCNCSSKWCRDHSRCC* 


TTX-r 


[55] 


PIMA 


ZRLCCGFOKSCRSRQCKOHRCC* 


Navl.4 


[38] 


GIIIA 


RDCCTOOKKCKDRQCK-OQRCCA* 


Navl.2, 1.4, 1.6-1.8 


[56] 


GIIIB 


RDCCTOORKCKDRRCK-OMKCCA* 


TTX-r 


[56] 


GIIIC 


RDCCTOOKKCKDRQCK-OLKCCA* 


TTX-r 


[56] 


CnlllA 


GRCCDVPNACSGRWCRDHAQCC* 


TTX-r 


[54] 


CnlllB 


ZGCCGEPNLCFTRWCRNNARCCRQQ 


TTX-r 


[54] 


SmIIIA 


ZRCCNGRRGCSSRWCRDHSRCC 


TTX-r, TTX-s 


[57] 


SIIIA 


ZNCCNGGCSSKWCRDHARCC* 


TTX-r, TTX-s 


[58] 


CIIIA 


GRCCEGPNGCSSRWCKDHARCC* 


TTX-r 


[59] 



doi:1 0.1 371 /journal.pone.00941 22.t001 



To confirm Tvl de novo assignment, the peptide was chemically 
synthesized using Fmoc Solid Phase Peptide Synthesis (SPPS). Due 
to inherent uncertainty in the amino acid sequence at position 3, 
both IleS and Leu3 containing peptides were synthesized and 
subjected to the same chemical derivatization strategy. The two 
synthesized Tvl teretoxins both resulted in an MS/MS spectrum 
identical to the native peptide, confirming the de novo sequence 
assignment in which TR is the order of the first two amino acids, 
but with remaining I/L uncertainty at position 3 due to the mass 
spectrometer used (Fig. 2 and Figure S2 in File SI). The novel 
peptide teretoxin, Tvl from venomous marine snail Terebra 
variegata was thus identified and its sequence was determined. 

Synthesis, Characterization and Oxidative Folding of Tvl 

Tvl was chemically synthesized using Fmoc Solid Phase Peptide 
Synthesis without prior knowledge of disulfide connectivity, and all 
cysteine residues were protected with a standard side chain 
protecting group, trityl. The synthesized Tvl peptide was purified 
to 99% using RP-HPLC (Figure S3 in File SI) and its mass was 
confirmed by MALDI mass spectrometry. The expected average 
mass of the linear peptide [Mh-H] 2316.3 m/z was observed 
(Figure S4 in File SI). Tvl was folded using thiol-assisted air 
oxidation reaction and purified by HPLC (Figure S3 in File SI). 
The mass of folded Tvl, 2310.7 m/z, is reduced by 6 Da 
confirming that three intramolecular disulfide bonds had been 
formed (Figure S5 in File SI). No differences were observed in 
synthesis and folding kinetics of the two Tvl peptides with He or 
Leu at the third amino acid position, and both versions of Tvl 
were used to determine the disulfide connectivity of the peptide. 



Disulfide Mapping of Tvl 

Due to the limited number of enzymatic cleavage sites and the 
presence of adjacent cysteine residues, traditional strategies for 
disulfide mapping such as tryptic digestion and acid hydrolysis 
were difficult to implement. The disuffide connectivity of Tv 1 was 
independently determined by MS /MS mapping and NMR 
spectroscopy. For MS/MS mapping, a partial reduction and dual 
alkylation protocol was applied using reducing agent TCEP-HCl 
(Tris(2carboxyethyl)phosphine hydrochloride) and alkylating 
agents NEM (N-ethylmaleimide) and lAM (iodoacetamide) [28- 
30]. Chemically synthesized and oxidatively folded Tvl peptide 
was first partially reduced at pH 3.0 with TCEP, and subsequently 
alkylated with NEM. The partially alkylated Tvl -NEM peptide 
was reduced further with TCEP at pH 7.0 and alkylated again 
with lAM this time. Dual NEM/IAM alkylation resulted in Tvl 
peptide species that were labeled with 2, 4 or 6 NEM and lAM 
groups (Fig. 3 and Figure S6 in File SI). The location of NEM and 
lAM modifications in each of the six partially reduced species was 
determined by matching the MS/MS b- and y- series ions to 
theoretical patterns. All six partially reduced species provided 
complementary, non-contradictory information and the resulting 
connectivity was determined to be Cys4— Cys20, Cys5-Cys21 and 
Cys7-Cysl6 (Fig. 3 and Table SI in File SI). Tvl has a cysteine 
scaffold similar to M Superfamily conotoxins (Table 1), but it 
displays a disulfide pattern previously unknown in native cone 
snail peptides [23,31-35] with a disulfide bridge formed between 
the two middle cysteines (Cys7 and Cysl6) and a parallel 2- 
disulfide bond "staple" (Cys4-Cys20 and Cys5-Cys21) linking the 
N- and C-termini of the peptide. 
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Figure 2. Sequence identification of Tvl teretoxin. Shown is the native toxin (blacl<) and the synthesized toxin (green) IVIS/MS spectrum 
recorded on a (M +6H)+6 ion after conversion of cysteine residues to dimethyl lysine analogs. Note the overlap of all major fragment ions. The 
sequence is given above the spectrum and observed c and z-type fragment ions are indicated in the sequence with 1 and [, respectively. Doubly 
charged fragment ions of type c and z are labeled with +2, triply charged ions are of type c and z are indicated with *, z-type fragment ions that 
resulted from cleavage at cysteine with subsequent loss of the cysteine side chain are denoted in italic [60] and charge reduced species are labeled in 
the spectrum with #. The spectrum was recorded at a resolution of 7500 at m/z 400 and all fragment ions have a mass accuracy of better than 
5 ppm. 

doi:1 0.1 371 /journal.pone.00941 22.g002 



Characterization of Tvl by NiVlR Solution Structure 
Derivation 

To confirm the disulfide bond connectivity derived from MS, 
the solution structure of Tvl was derived using standard 
homonuclear proton NMR techniques on unlabeled folded 
synthetic peptide. Initial NMR experiments were carried out on 
Tv 1 peptide containing either He or Leu at the third position. The 
spectra were very similar for both versions of the peptide, with the 
only differences found in the signals of the Ile3 and Leu3 residues 
themselves, indicating that both peptides adopt the same structure. 
Subsequent NMR structural calculations were carried out using 
only the Ile3 form of the peptide. 

Proton assignments were obtained from 2D NOESY and 
TOCSY spectra, and carbon chemical shifts were assigned with 
the help of a natural-abundance C-HSQC spectrum (Table S2 
in File SI). Initial structure calculation runs were carried out using 
only manually assigned Nuclear Overhauser Effect (NOE) distance 
restraints. Disulfide connectivities were then determined based on 
proximity of cysteine residues in the 1 0 lowest-energy structures 
and were in agreement with the disulfide bond pattern derived by 
MS analysis. 

Direct evidence for the Cys7-Cysl6 disulfide bond was also 
observed from NOE crosspeaks between the H'^ of Cys 1 6 and the 
of Cys7, as well as a number of connectivities to the two 
flanking residues, Gly6 and Tyr8 (Figure S7 in File SI). Once the 
disulfide linkages were determined, additional distance restraints 
for the disulfide bonds were introduced. The final round of 
structure calculations was run with covalent bonds between the 
connected cysteine side chains. A bundle of the 10 lowest-energy 
structures was obtained with an RMSD of 0.42 A for backbone 



and 0.74 A for all heavy atoms (Figure S8 and Table S3 in File 
SI). A query for other peptide structures similar to the determined 
Tvl structure using the software packages VAST [36] or 
PDBeFOLD [37] did not yield any hits, indicating that the 
structure of Tvl represents a new fold for small peptide toxins. 
The most prominent and best-defined structural feature is a (3- 
hairpin from Cys 7 to Cys 16 that is clamped together by the 
disulfide bond formed between these two residues (Fig. 1 and 
Figure S9 in File SI). The remainder of the peptide wraps around 
the side chain aromatic ring of Tyr8, in particular through the 
formation of hydrophobic interactions with the side chains of Ile3 
and Vail 5. The N- and C-terminal loops are clamped together in 
an antiparallel way through the Cys4-Cys20, Cys5-Cys2 1 double- 
disulfide bond arrangement, giving the whole teretoxin the shape 
of a flattened ellipsoid (Figure S8 in File SI). Coordinates for the 
Tvl structures have been submitted to the PDB and are available 
under accession code 2mrx. 

Although Tvl has a cysteine scaffold similar to M Superfanuly 
conotoxiiis (Table 1), the observed differences in fold and disulfide 
bonding pattern are not surprising given the lack of sequence 
homology to M superfamily conotoxins. It is interesting to note 
that the loop between C3 and C4 is at least 4 amino acid residues 
longer in Tvl compared to any of the other listed conotoxins, 
which is in agreement with the fact that these are the residues 
involved in the formation of the pi-hairpin that is lacking in 
conotoxins. No fi-strand structures are fovmd in any of the other 
M-superfamily conotoxins (Fig. 4) [28-32]. To our knowledge, the 
only other example of Tvl's disulfide bond pattern was found in 
an active non-native minor refolding product for )i-conotoxin 
PIIIA [38] . The 3D structures of both the native and non-natively 
refolded forms of PIIIA are very different from the structure of 
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Figure 3. Disulfide mapping of Tvl. a) RP-HPLC fractionation of all differentially folded peaks (on HPLC gradient of 20% to 35% buffer B (80% 
Acetonitrile, 0.1% TFA) in buffer A (water, 0.1% TFA) over 75 min. b) MS-IVIS spectrum of partially reduced peptide peaks showing representative 
alkylated b and y ions, c) Schematic representation of Tvl with predicted disulfide bridges. 
doi:1 0.1 371 /journal.pone.00941 22.g003 



Tvl as they contain segments of a-helix and no fi-sheet structure. 
The difiFerence between Tvl and the PIIIA refolded product with 
the same disulfide bridge arrangement further supports the 



conclusion that Tvl 
of peptide toxins. 



representative of a new and unique group 






Tvl 



SmIIIA 



mr3e 



KIIIA 



Figure 4. Structural comparison of Tvl teretoxin and M-superfamily conotoxins. Comparison of the NMR structure of Tvl with that of M 
superfamily conotoxins SmIIIA, mr3e, and KIIIA reveals significant structural differences between Tvl and these conotoxins despite all having the 
same cysteine scaffold CC-C-C-CC. All structures are shown in cartoon representation with disulfide bonds highlighted in yellow. All figures were 
prepared using PyMol (www.pymol.org). Conotoxin structural references are as follows: Mrllle [33] SmIIIA [61] and KIIIA [62]. 
doi:1 0.1 371/journal.pone.00941 22.g004 
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Bioactivity Assay 

To determine if synthetic Tvl is biologically active, the peptide 
was injected into polychaete worms, which are the natural prey of 
terebrid snails. Injections of 20 |iM Tvl (He) into N. virem 
polychaete worms caused partial paralysis (Fig. If, Table 2). This 
fmding is consistent with the recent study where crude terebrid 
venom extract exhibited nAChR inhibitory activity [39]. Inhibi- 
tion of nicotinic receptors at the neuromuscular junction can result 
in paralysis. Further experiments are underway to determine the 
specific molecular target of Tvl. 

Conclusion 

A novel teretoxin, Tvl, from Terebra variegata was structurally 
identified and found to be biologically active. The strategy applied 
to the discovery and characterization of Tvl allows for the 
unequivocal analysis of minuscule amounts of sample. Four 
nanograms of venom sample was used to elucidate the amino 
acid sequence of Tvl. Such a minute amount of sample would 
previously be inaccessible for natural product discovery from 
venomous organisms. There are over 400 different species of 
terebrids, each predicted to express between 50-200 difiFerent 
peptide toxins in their venom. With an integrated approach of 
mass spectrometry, chemical peptide synthesis, and NMR as 
applied to Tvl, it will be possible to access the substantial quantity 
of disulfide-rich peptides (>80,000) in the Terebridae. The novel 
structural features of Tvl suggest that teretoxins will reveal 
different mechanisms of action and different target specificities 
compared to disulfide-rich peptide toxins previously characterized 
and used as tools in neuroscience and as drug leads. In addition, 
an abundance of gastropods, including most species of the 
venomous marine conoidean snails, are smaller than 8 mm in 
shell size [40], highlighting the need for an efficient microscale 
mass spectrometry sequencing strategy. The sample limited 
integrated approach presented here can significantly enhance 
discovery and characterization of the vast disulfide-rich natural 
product peptide libraries produced by various venomous organ- 
isms. 

Methods 

Venom Extraction and De-Novo Sequencing 

Four Terebra variegata snails were collected by MH from the Las 
Perlas Archipelago in Panama on an expedition with the RV 
URRACA and glands were dissected. Scientific research permits 
were provided by the Smithsonian Tropical Research Institute 
Permit Office (STRI-SPO) and The Panama Aquatic Resources 
Authority (ARAP). Pooled crude venom contents from all 4 
specimens were extracted with aid of light sonication in 30% 
acetonitrUe/water acidified with 0.1% trifluoroacetic acid and 
centrifuged. The supernatant was lyophilized and stored at — 

Table 2. Tvl bioactivity in polychaete wornns. 



20°C. An aliquot of crude venom was pressure loaded onto a self- 
packed PicoFrit column (New Objective) with integrated emitter 
tip (360 |im o.d., 75 |im i.d., 15 |im tip), packed with 6 cm of 
reverse-phase C18 material (Alltima C18 5 |.im beads, AUtech 
Associates), rinsed for 10 min with 0.1 M acetic acid and 
subsequently gradient eluted with a linear gradient from 0 to 
100% B in 30 min (A = 0.1 M acetic acid, B = 70% acetonitrUe in 
0.1 M acetic acid, flow rate 120 nL/min) into an LTQ^Orbitrap 
XL with ETD (ThermoFisher Scientific), using a home-built micro 
electrospray source with a liquid junction. 

The instrument was operated in data dependent mode cycling 
through fuU scan (300-2,000 m/z, single scan) followed by 8 MS/ 
MS scans on the 8 most abundant ions from the immediate 
preceding full scan. Next, aliquots of reduced and iodoacetamide 
alkylated crude venom were pressure loaded onto a self-packed 
PicoFrit column (New Objective) with integrated emitter tip as 
described above and eluted using the same gradient as above into 
the LTQ_ Orbitrap XL with ETD. Instrument was operated in a 
data dependent mode cycling through a fuU scan (300-2,000 m/z, 
single scan) foUowed by 8 MS/MS scans on the 4 most abundant 
ions from immediate preceding fuU scan. The instrument was 
programmed to acquire first 4 CAD and then 4 ETD spectra on 
the 4 most abundant ions. The cations were isolated with a 2-Da 
mass window and set on a dynamic exclusion list for 30 seconds 
after they were first selected for MS/MS. Target value for MS/ 
MS was set to 10* cations and 10^ anions of fluoranthene. For 
CAD an activation time of 30 ms was chosen. For ETD the ion/ 
ion reaction time was set to 100 ms. In both cases 4 scans were 
averaged per MS/MS spectrum. In a majority of cases, complete 
sequence elucidation after this initial analysis was not possible, but 
based on this initial screen optimal charge states for subsequent 
targeted ETD analysis were chosen. For subsequent analysis either 
reduced, iodoacteamide alkylated, or dimethyl lysine analog 
versions of toxin were pressure loaded onto the self-packed 
columns as described above. The column was rinsed for 10 min 
(for reduced or iodoacetamide alkylated toxins) or for 20 min (for 
lysine dimethyl analog versions of toxins) with 0. 1 M acetic acid 
and gradient eluted into the LTQ^ XL as described above. To 
obtain higher quality MS/MS spectra in this second analysis, 12 
scans were averaged and the data acquired in high resolution with 
7,500 resolution setting. The spectra were manually interpreted 
(Figures SI and S2 in File SI). 

Synthesis and Purification of Tvl 

Tvl peptide was synthesized by microwave assisted Fmoc SPPS 
on a CEM Liberty synthesizer using standard side chain 
protection. Standard Fmoc cysteine with side chain protecting 
group Trityl (Fmoc-Cys (Trt)-OH) was used for all cysteine 
residues in the Tvl peptide. FoUowing treatment of peptidyl resin 
with Reagent K [92.5% TEA (Trifluoroacetic acid), 2.5% TIS 



Activity Time spent 

Tvl (20 mM) NSS 

atO 0.995 1.45 

Between 0 and 1 74.47 0 

Between 1 and 2 19.98 0 

Between 2 and 3 4.55 98.55 

doi:l 0.1 371 /journal.pone.00941 22.t002 
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(Triisopropylsilane), 2.5% EDT (1,2 Ethanedithiol) and 2.5% 
water, 4 hr)] and cold ether precipitation, crude Tvl was purified 
by RP-HPLC using an X-Bridge semipreparative column (10 X 150 
mm, 5 [im particle size. Waters Corporation, MUford, MA, USA). 
Elution was carried out at 5 mL/min with 20% bufTer B (80% 
Acetonitrile, 0.1% TFA) and 80% buffer A (0.1% TEA) for the 
first 5 min, then increasing buffer B to 35% in 45 min. HPLC 
buffer composition remained the same all through the work. Purity 
of collected peptide was confirmed by RP-UHPLC using an 
Acquity UPLC (BEH 300 CIS 1.7 |xm. Waters Corporation) 
column and eluted using a linear gradient from 0% to 75"/(> buff-r 
B in 3.5 min (Figure S3 in File SI). The identity of .synthesized 
peptide was confirmed by molecular mass measurement of 
purified peptide using MALDI-TOF (Waters, Micromass, 
CHaCN matrix) (Figure S4 in File SI). 

Oxidative Folding of Tvl 

A one-step thiol-assisted oxidation was used to prepare folded 
Tvl peptide. The linear peptide (20 (tM) was incubated in O.I M 
Tris -HCl, 0.1 M NaCl, 100 nM EDT A, 1 mM GSH, 1 mM 
GSSG, pH 7.5. The folding reaction was terminated by acidifi- 
cation with 8% formic acid at 15 min, 30 min, 1, 2, 3, 4, and 24 h 
and the folding yield monit()r(;d using UHPLC. A preparative 
scale folding reaction was then conducted at an optimized time of 
2 h, and the folded peptide was purified using X-Bridge 
semipreparative column (as mentioned earlier). Elution was 
carried out at 5 mL/min with 15% buffer B and 85% buffer A 
for the first 5 min, then increasing buffer B to 35% in 45 min. The 
purity was confirmed using Acquity UHPLC column (Figure S3 in 
File SI), and the molecular mass of the oxidized peptide was 
confirmed by MALDTTOF (Figure S5 in File SI). 

Disulfide Mapping by Tandem MS/MS 

A partial reduction and dual alkylation strategy was employed 
using reducing agent TCEP (Tris carboxy ethyl-phosphine 
hydrochloride, Thermo-Fisher) and alkylating agents NEM (N- 
ethyknaleimide, Sigma) and lAM (iodoacetamide, Sigma). Con- 
ditions for the partial reduction and alkylation of folded peptides 
were optimized modifying previous protocols [30,41]. Briefly, 9 |il 
of folded peptide (1 mM) was dissolved in 0.1 M citrate buffer, 
pH 3.0 and 1 Hi of TCEP (20 mM in 0. 1 M citrate buffer, pH 3.0) 
was added to it. Reduction was carried out at 4°C for 90 min and 
was terminated by injecting onto an RP-HPLC column [X-Bridge 
analytical column (5x150 mm, 5 |Xm particle size. Waters 
Corporation, Milford, AL\, USA)] and eluted with 2(1% buffer B 
in buffer A for 5 min, then increasing buffer B to 35% in 75 min. 
Six partially reduced peaks (Figure S6 in File SI) were each 
collected direcdy into microcentrifuge tubes containing 100 |4l 
NEM (200 mM in 0.1 M citrate buffer pH 3) and incubated at 
37°C for 2 h in the dark. Partially reduced and alkylated reactions 
were then lyophilized and desalted to remove unreacted NEM 
using UHPLC. The desalted fractions were incubated with 
100 mM TCEP at 55°C for 1 h to completely reduce the NEM 
labeled peptide and alkylated with 275 nM lAM for 2 h at 37°C 
in dark. The alkylation steps resulted in peptide species that were 
labeled with 2, 4 or 6 NEM or 1AM groups. The labeled 
ly()pliiliz(;d species were suspended in 0.5% formic acid and 
subjected to LC-MS/MS analysis on 43 mm HPLC-Chip/Q; 
TOF using a 7 min gradient to resolve fully labeled species from 
those with incomplete alkylations. Data was collected in targeted 
and auto MS/MS mode and processed using molecular feature 
extraction (MFE) software to detect unique peptide features (MS 
and MS/MS spectra) followed by sequence matching on 
MassHunter Bioconfirm Qual B.05 software (Table SI in File SI). 



NMR Solution Structure Derivation 

Samples for Nuclear Magnetic Resonance (NMR) studies were 
prepared by dissolving lyophilized oxidized peptide into either a 
9:1 ratio of H20:D20 or in 100% DjO at a concentration of 
—260 |xM. Sample pH was adjusted to pH 6.0 and 50 ^M 
trimethylsilyl propanoic acid was added as a reference compound. 
Spectral acquisition was carried out at 5°C on a Varian Inova 
600 MHz NMR spectrometer equipped with a cryogenically 
cooled HCN probe. 2D-TOCSY spectra were acquired in 90"/o 
H2O (64 scans with 1024x512 points) and D2O (48 scans with 
1024x512 points) with spin-lock times of 60 ms, whereas 2D- 
NOESY spectra were acquired with mixing times of 200 ms in 
HjO with a matrix size of 1024x512 points in 64 scans and in 
D2O with a similar matrix size in 1 12 scans. A natural abundance 
'h'^c hsqc correlation spectrum was also acquired in 1024 
scans with a matrix size of 1024x80 points. Acquired spectra were 
processed with NMRPipe [42] and analyzed with CCPNmr 
Analysis [43]. Proton assignments were made by comparison of 
TOCSY and NOESY spectra, and used to assign the ''^C chemical 
shifts through the HSQC spectrum (Table S2 in File SI). Distance 
restraints for structural calculations were obtained from 2 NOESY 
spectra: one spectrum collected in 90% H2O with a mixing time of 
200 ms and another spectrum collected in D2O with a mixing 
time of 400 ms. Manual assignment was performed for about 90% 
of all NOE cross-peaks. Dihedral angle restraints were derived 
from the assigned chemical shifts using the software TALOS+ 
[44]. Structural calculations were carried out with ARLA2 [45]/ 
CNS [46] . Initial structural calculations were conducted using only 
the manual NOE assignments (with distances automatically 
calibrated by the software based on peak volumes) and without 
any restraints for the disulfide bonds. An ARIA2 run consisting of 
8 iterations of simulated annealing followed by explicit water 
refinement of the 10 lowest-energy structures was performed. 
Disulfide connectivities were determined based on this initial 
structural bundle and included in the final round of structural 
calculations as covaleiit bonds. Eight itc'rations were performed in 
ARIA with 20 structures generated per iteration and the 7 lowest 
energy structures of each iteration forming the structural ensemble 
for the next iteration. For this run, NOE assignments were allowed 
to be automatically adjusted by the software during the 
calculation. For the last iteration, 100 structures were generated, 
and ten lowest-energy solutions were subjected to explicit solvent 
refinement. The quality of the final structural bundle was assessed 
by Procheck NMR (Table S3 in FUe SI) [47]. 

Bioactivity Assay on Polychaete Worms 

Mreis virens (polychaete worms) (<3.5 g) were maintained at 4°C 
in salt water prior to injection. Worms were laid on a Styrofoam 
block and injected with an insulin syringe. A modified needle cap 
was used to maintain a consistent 1.5 mm depth of needle 
puncture. Worms were injected between the 5* and 7* segments 
on the ventral anterior end with the aim of targeting the ventral 
nerve cord. Control worms were injected with 12 ^1 of 0.9% 
normal saline solution (NSS). Lyopholized samples of Tvl toxin 
were resuspended in NSS and diluted to a concentration of 20 hM 
which was confirmed using a ThermoScientific Nanodrop 2000c 
spectrophotometer. At least three worms were used for each trial- 
one saline control and two injected with Tvl toxin for a total of 3 
trials. After injection, the worms were transferred to side-by-side 
containers filled with room temperature salt water. The worm 
movements were then recorded by video in order to evaluate the 
effects of the toxin over a period of at least 3.5 hours. Varying 
temperature ehcits different phenotypic responses from N. virens: 
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water at 4°C has an anesthetizing efiect, while room temperature 
prompts a vigorous response [48] . 

Taking the worms from 4°C water and immersing them in 
room temperature water (25°C) immediately after injection 
stimulates an active response from the worms, making it easier 
to recognize paralytic (or partially paralytic) effects induced by a 
novel toxin in the first minutes after injection. Worms were 
observed every other day for two weeks after injection. The 
recorded videos were analyzed blind by two independent 
observers; the activity of each worm over a small time interval 
was noted on a scale of 1-3, with 3 being extremely active and 1 
being inactive. Each video was analyzed and annotated three times 
by each observer. Observations for each worm began at the time 
of injection (Table S4 in File SI). The numerical results of the two 
independent viewings were averaged together if results were 
similar, if two viewers disagreed, a third obser\'er was employed to 
review the disputed section of video. Results of relative activit)' vs. 
time were plotted to quantify the phenotypic eflects of the toxin 
and saline injection on jV. virens (Fig. 11). The data was also used to 
calculate what percentage of time each worm spent inactive (a 
rating of 1), active (an averaged rating of 1.5—2), and extremely 
active (rating 2.5-3) (Table 2). 

Supporting Information 

File SI This file contains Figure SI— Figure S9 and Table 
Si-Table S4. Figure SI, CAD of native Tvl. MS/MS spectrum 
recorded on a (M +2H)"^^ ion after reduction of cysteine residues. 
The sequence is given above the spectrum and observed b, a and 
y-type fragment ions are labeled in the spectrum. Observed 
peptide backbone cleavage is indicated in the sequence above with 
1 and |_ for N- and C-terminal fragment ions, respectively. Doubly 
charged fragment ions are labeled with The neutral loss of 
water from the precursor ion is shown as [M+2H]"^2-H20, but 
neutral losses of fragment ions are not labeled. The spectrum was 
recorded at a resolution of 7500 at m/ z 400 and all fragment ions 
have a mass accuracy of better than 5 ppm. Figure S2, ETD of 
native (black) and synthetic Tvl (blue). MS/MS spectrum 
recorded on a (M +6H)"'"^ ion after conversion of cysteine residues 
to dimethyl lysine analogs. The sequence is given above the 
spectrum and observed c and z-type fragment ions are indicated in 
the sequence with "| and [ , respectively. Doubly charged fragment 
ions of type c and z- are labeled with triply charged ions are of 
type c and z are indicated with *, z-type fragment ions that 
resulted from cleavage at cysteine with subsequent loss of the 
cysteine side chain are denoted in italic and charge reduced species 
are labeled in the spectrum with The spectrum was recorded at 
a resolution of 7500 at m/z 400 and all fragment ions have a mass 

References 

1. Koh CY, Kini RM (2012) From snaJce venom toxins to therapeutics- 
cardiovascular examples. Toxicon 59: 497-506. 

2. King GF, Hardy MC (2013) Spider-venom peptides: structure, pharmacology, 
and potential for control of insect pests. Annu Rev Entomol 58: 475—496. 

3. Vetter I, Lewis RJ (2012) Therapeutic Potential of Cone Snail Venom Peptides 
(Conopeptides). Curr Top Med Chem 12: 1546-1552. 

4. King Gf (201 1) Venoms as a platform for human drugs: translating toxins into 
therapeutics. Expert Opin Biol Ther 11: 1469-1484. 

5. Cheek S, Krishna SS, Grishin NV (2006) Structural classification of small, 
disulfide-rich protein domains. J Mol Biol 359: 215-237. 

6. Cushman DW, Ondetti MA (1991) History of the design of captopril and related 
inhibitors of angiotensin converting enzyme. Hypertension 17: 589-592. 

7. Miljanich GP (2004) Ziconotide: Neuronal Calcium Chaimel Blocker for 
Treating Severe Chronic Pain. Curr Med Chem 11: 3029-3040. 

8. Prashanth JR, Lewis RJ, Dutcrtre S (2012) Towards an integrated venomics 
approach lor accelerated eonopeptide discovery. Toxicon 60: 470-477. 

9. Dutertre S, Lewis RJ (2010) Use of venom peptides to probe ion channel 
structure and function. J Biol Chem 285: 13315-13320. 



accuracy of better than 5 ppm. Figure S3, RP-UHPLC chro- 
matograms of Tvl hnear and oxidized peptide at 214nm. During a 
pilot folding reaction, over 90% of the linear Tvl peptide fully 
oxidized and showed a peak at 1.58 minute in comparison to 1.83 
minute linear Tvl peak at the gradient of 0-75'% buffer B (80% 
acetonitrUe, 0. 1 % TFA) in buffer A (0. 1 % TFA) within two hours. 
Figure S4, Analysis of linear Tvl peptide by MALDI-TOF mass 
spectrometry. MALDI-TOF spectrum of Tvl peptide using a- 
Cy2Lno-4~hydrox)>cinnamic acid matrix. Figure S5, Analysis of 
oxidized Tvl peptide by MALDI-TOF mass spectrometry. 
MALDI-TOF spectrum of Tvl peptide using a-Cyano-4-A)'(/ro- 
xycinnamic acid vaMxhL. Figure S6, RP-UHPLC analysis of partially 
reduced Tvl peptide. On a UPLC (BEH 300 CI 8 1.7 jim. Waters 
Corporation, Mjlford, MA, USA) column, a linear gradient of 0— 
75% buffer B (80% acetonitrile, 0.1% TFA) in buffer A (0.1% 
TFA) over 6 minutes and peaks were assigned by their degree of 
reduction. Labels indicate how many disulfides are present. Figure 

57, NOE contacts confirming the C7-C16 disulfide bond. An 
overlay of the H'^H" fingerprint region of the NOESY (in black) 
and TOCSY (in blue) spectra shows NOE crosspeaks linking Cys 
16 and Cys 7 as well as contacts in the residues flanking the C7- 
C16 disulfide bond (C7-S17, Y8-C16, G6-C16, G6-N18). Figure 

58, Bundle of the 10 lowest energy- structures of Tvl after explicit 
water refinement. The 10 lowest energy structures are shown in 
stick representation, displaying the tight convergence found in the 
final structural bundle. Figure S9, Structure of Tvl. An overlay of 
a cartoon representation and a stick model of the lowest-energy 
structure of Tvl shows the P-sheet character of the peptide and 
reveals the important role of the Tyr 8 side chain in the formation 
of a small hydrophobic core. Table SI, Predicted and observed of 
b and y ions of differentially alkylated peptides by auto and 
targeted MS/MS analysis. Table S2, Chemical shift assignments 
of Tvl (in ppm). Table S3, Structural statistics for the final 10 
models of Tvl. Table S4, Sample raw data of Tvl bioactivity in 
polychaete worms. 

(DOC) 

Acltnowledgments 

The authors thank Dr. Vadiraj Bhat, Bryan Lavery, and Doug Postl from 
AgUcnt Technologies, and Dr. Jianqin Zhuang, NMR facility manager of 
the College of Staten Island, for help with running experiments. 

Author Contributions 

Conceived and designed the experiments: PA MH. Performed the 
experiments: PA AG MB BU VR JQ. PM. Analyzed the data: PA AG 
MB BU JQ. PM BC SP. Contributed reagents/materials/analysis tools: 
MH BU SP BC. Wrote the paper: PA MH BU SP BC. 



10. Favreau P, Laure Menin SM, Perret F, Cheneval O, StOcUin M, et al. (2006) 
Mass spectrometry strategies for venom mapping and peptide sequencing from 
crude venoms: Case applications with single arthropod specimen. Toxicon 47: 
676-687. 

11. Favreau P, Stocklin R (2009) Marine snail venoms; use and trends in receptor 
and channel neuropharmacology. Curr Opinion Pharmacol 9: 594—601. 

12. Escoubas P, Qtiinton L, Nicholson (iM (2008) Venomics: unravelling the 
complexity' of animal \'enoms with mass spectrometry. J Mass Spectrom 43: 
279-295. 

13. Vetter 1, Davis JL, Rash LD, .\nangi R, Mobli M, et al. (201 1) Venomics: a new 
paradigm for natural products-based drug discovery. Amino Acids 40: 15—28. 

14. Quinton L, Demeure K, Dobson R, Gilles N, De Pauw E (2007) New Method 
for Characterizing Highly Disulfide-Bridged Peptides in Complex Mixtures: 
Application to Toxin Identification from Crude research articles. J Proteome 
Res 6: 3216-3223. 

15. Tayo LL, Lu B, Cruz LJ, Yates JR III (2010) Proteomic analysis provides insights 
on venom processing in Conus textile. J Proteome Res 9: 2292-2301. 



PLOS ONE I www.plosone.org 



8 



April 2014 | Volume 9 | Issue 4 | e94122 



Sample Limited Venom Peptide Cliaracterization 



16. Bhattachary\'-a M, Gupta K, (iowd KH, Balaram P (2013) Rapid mass 
spectromctric determination of disuU'idc connectivity in peptides and proteins. 
Mol Biosyst 9: 1340-1350. 

17. Bhatia S, Kil YJ, Ueberheide B, Chait BT, Tayo L, et al. (2012) Constrained de 
novo sequencing of conotoxins. J Proteome Res 11: 4191^200. 

18. Castelin M, Puillandre N, Kantor YI, Modica MV, Terryn Y, et al. (2012) 
Macroevolution of venom apparatus innovations in auger snails (Gastropoda; 
(.^onoidea; Terebridac). Mol Phylogenet Evol 64: 21-44. 

19. Holford M, Puillandre N, lcrr\'-n Y, Cruaud G, Olivera B, et al. (2009) 
Evolution of the Toxoglossa venom apparatus as inferred by molecular 
phylogeny of the Terebridac. Mol Biol Evol 26: 15—25. 

20. Olivera BM, Teichert RW (2007) Diversity of the neurotoxic Conus peptides. 
Mol Interv 7: 251-260. 

21. Olivera BM (2006) Gonus peptides: biodiversity-based discover)' and exoge- 
nomics.J Biol Ghem 281: 31173-31177. 

22. Terlau H, Olivera BM (2004) Gonus venoms: a rich source of novel ion channel- 
targeted peptides. Physiol Rev 84: 41-68. 

23. fian Y-H, Wang Jiang H, Liu L, Xiao C, et al. (2006) Characterization of 
novel M-superfamily conotoxins with new disulfide linkage. FEBS J 273: 4972— 
4982. 

24. Bulaj G (2008) Integrating the discovery pipeline for novel compounds targeting 
ion channels. Curr Opin Chem Biol 12: 441^47. 

25. Imperial JS, Watkins M, Chen P, Hillyard DR, Cruz LJ, et al. (2003) The 
augertoxins: biochemical characterization of venom components from the 
toxoglossate gastropod Terebra subulnta. Toxieon 41: 391-398. 

26. Imperial JS, Kantor Y, Watkins M, lii EMH, Stevenson B, et al. (2007) 
Venomous Auger Snail Hastula (Impagcs ) hcctica (Linnaeus, 1758 ): Molecular 
Phylogeny, Foregut Anatomy and Comparative Toxinology. 308B: 744—756. 

27. Ueberheide BM, Fenyd D, Alewood PF, Chait BT (2009) Rapid sensitive 
analysis of cysteine rich peptide venom components. Proc Natl Acad Sci USA 
106: 6910-6915. 

28. Nair SS, Nilsson GL, Emmett MR, Schaub TM, Gowd KH, et al. (2006) De 
Novo Sequencing and Disulfide Mapping of a Bromotryptophan-Containing 
Gonotoxin by Fourier Transform Ion Gyclotron Resonance Mass Spectrometry. 
Anal Ghem 78: 8082-8088. 

29. Jakubowski JA, Sweedler JV (2004) Sequencing and mass profiling highly 
modified conotoxins using global reduction/ alkylation followed by mass 
spectrometry. Anal Chem 76: 6541-6547. 

30. Goransson U, Craik DJ (2003) Disulfide Mapping of the Cyclotide Kalata Bl. 
J Biol Chem 278: 48188-48196. 

31. Corpuz GP, Jacobscn RB, Jimenez EC, Watkins M, Walker G, et al. (2005) 
Definition of the M-conotoxin superfamily: characterization of novel peptides 
from moUuseivorous Gonus venoms. Biochemistry 44: 8176-8186. 

32. Santos AD, Mcintosh JM, Hillyard DR, CJruz LJ, Olivera BM (2004) The A- 
superfamily of conotoxins: structural and functional divergence. J Biol Chem 
279: 17596-17606. 

33. Du W-H, Han Y-H, Huang F-J, Li J, Chi C-W, et al. (2007) Solution structure of 
an M-1 conotoxin with a novel disulfide linkage. FEBS J 274: 2596-2602. 

34. McDougal CM, Poulter CD (2004) Three-dimensional structure of the mini-M 
conotoxin mr3a. Biochemistry 43: 425—429. 

35. Jacob RB, McDougal OM (2010) The M-superfamily of conotoxins: a review. 
Cell Mol Life Sci 67: 17-27. 

36. Gibrat J-F, Madej T, Bryant SH (1996) Surprising similarities in structure 
comparison. Curr Opin Struct Biol 6: 377-385. 

37. Krissinel E, Henrick K (2004) Secondary-structure matching (SSM), a new tool 
for fast protein structure alignment in three dimensions. Acta Crystallogr D Biol 
Crystallogr 60: 2256-2268. 

38. Tietze AA, Tietze D, Ohlenschlager O, Leipold E, Ullrich F, et al. (2012) 
Structurally diverse |l-conotoxin PIIIA isomers block sodium channel NaV 1.4. 
Angew Ghem Int Ed Engl 51: 4058-4061. 

39. Kendel Y, Melaun G, Kurz A, Nicke A, Peigneur S, et al. (2013) Venomous 
secretions from marine snails of the Terebridae family target acetylcholine 
receptors. Toxins (Basel) 5: 1043-1050. 

40. Bouchet P, Lozouet P, Sysoev AV (2009) An inordinate fondness for turrids. 
Deep-Sea Res. II 56: 1724-1731. 



41. Foley SF, Sun Y, Zheng TS, Wen D (2008) Picomole-level mapping of protein 
disulfides by mass spectrometry following partial reduction and alkylation. Anal 
Biochem 377: 95-104. 

42. DelagUo F, Grzesiek S, Vuister GW, Zhu G, Pfeifer J, et al. (1995) NMRPipe: a 
multidimensional spectrjil processing system based on UNIX pipes. J Biomol 
NMR 6: 277-293. 

43. Vranken WF, Boucher W, Stevens TJ, Fogh RH, Pajon A, et al. (2005) The 
GGPN data model for NMR spectroscopy: development of a software pipehne. 

Proteins 59: 687-696. 

44. Shen Y, Delaglio F, (..ornilescu G, Bax A (2009) TALOS-h: a hybrid method for 
predicting protein backbone torsion angles from NMR chemical shifts. J Biomol 
NMR 44: 213-223. 

45. Rieping W, Habeck M, Bardiaux B, Bernard A, Malliavin TE, et al. (2007) 
ARIA2: automated NOE assigmnent and data integration in NMR structure 
calculation. Bioinformatics 23: 381—382. 

46. Brtinger AT, Adams PD, Clore GM, DeLano WL, Gros P, et al. (1998) 
Crystallography & NMR System: A New Software Suite for Macromolecular 
Structure Determination. Acta Grystallogr D Biol (.:r\"staIlogr D54: 905-921. 

47. Laskowski RA, Rullmann JAC, MacAnhiir M\\\ Kaptcin R. Thornton JM 
(1996) AQUA and PROCJHEGK-NMR: Programs for checking the quality of 
protein structures solved by NMR. J Biomol NMR 8: 477-486. 

48. ImperialJS, Chen P, Sporning A, Terlau H, Daly NL, et al. (2008) Tyrosine-rich 
conopeptides affect voltage-gated K-h channels. J Biol Ghem 283: 23026—23032. 

49. Wang Q; Jiang H, Han Y-H, Yuan D-D, Chi G-W (2008) Two different groups 
of signal sequence in M-superfamily conotoxins. Toxieon 51: 813-822. 

50. Kaas Q, Westermann J-C, Halai R, Wang CKL, Craik DJ (2008) ConoServer, a 
database for conopeptide sequences and structures. Bioinformatics 24: 445—446. 

51. MeDougal OM, Turner MW, Ormond Aj, Poulter CD (2008) Three- 
dimensional structure of conotoxin tx3a: An m-1 branch peptide of the M- 
superfamily. Biochemistry' 47: 2826-2832. 

52. Franco A, Pisarewicz K, MoUcr G, Mora D, Fields GB, et al. (2006) 
Hyperhydroxylation: anew strategy for neuronal targeting by venomous marine 
molluscs. In MuUer E.G. Werner, editor. Progress in Molecular and Subcellular 
Biology. Springer, pp. 83-99. 

53. Holford M, Zhang M-M, Gowd KH, Azam L, Green BR, et al. (2009) Pruning 
nature: Biodiversity-derived discovery of novel sodium channel blocking 
conotoxins from Conm bullatus. Toxieon 53: 90-98. 

54. Zhang M-M, Fiedler B, Green BR, Gatlin P, Watkins M, et al. (2006) Structural 
and functional diversifies among mu-conotoxins targeting TTX-resistant sodium 
channels. Bioehemistrv 45: 3723-3732. 

55. Zhang M-M, Green BR, Gatlin P, Fiedler B, Azam L, et al. (2007) Structure/ 
function characterization of micro-conotoxin KIIIA, an analgesic, nearly 
irreversible blocker of mEunmahan neuronjil sodium channels. J Biol Chem 
282: 30699-30706. 

56. Adams DJ, Alewood PF, Craik DJ, Drinkwater RD, Lewis RJ (1999) Conotoxins 
and their potential pharmaceutical applications. Drug Dev Res 46: 219-234. 

57. West PJ, Bulaj G, Garrett JE, Olivera BM, Yoshikami D (2002) Mu-conotoxin 
SmlllA, a potent inhibitor of tetrodotoxin- resistant sodium channels in 
amphibian sympathetic and sensory neurons. Biochemistry 41: 15388-15393. 

58. Yao S, Zhang M-M, Yoshikami D, Azam L, Olivera BM, et al. (2008) Structure, 
dynamics, and selectivity of the sodium channel blocker mu-conotoxin SlllA. 
Biochemistry 47: 10940-10949. 

59. Shon KJ, Olivera BM, Watkins M, Jacobsen RB, Gray WR, et al. (1998) mu- 
Conotoxin PIIIA, a new peptide for discriminating among tetrodotoxin-sensitive 
Na channel subtypes. J Neurosci 18: 4473-4481. 

60. Ghalkley RJ, Brinkworth CS, Burlingame AL (2006) Side-chain fragmentation of 
alkylated cysteine residues in electron capture dissociation mass spectrometry. 
J Am Soc Mass Spectrom 17: 1271-1274. 

61. Keizer DW, West PJ, Lee EF, Yoshikami D, Olivera BM, et al. (2003) Structural 
basis for tetrodotoxin-resistant sodium channel binding by mu-conotoxin 
SmlllA. J Biol Ghem 278: 46805-46813. 

62. Khoo KK, Gupta K, Green BR, Zhang M-M, Watkins M, et al. (2012) Distinct 
disulfide isomers of n-conotoxins KIIIA and KIIIB block voltage-gated sodium 
channels. Biochemistry 51: 9826-9835. 



PLCS ONE I www.plosone.org 



9 



April 2014 | Volume 9 | Issue 4 | e94122 



